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Characterizing the Larkin-Ovchinnikov-Fulde-Ferrel phase induced by the 
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We discuss possible destinations from the chromomagnetic instability in color superconductors 
with Fermi surface mismatch Sfi. In the two-flavor superconducting (2SC) phase we calculate the 
effective potential for color vector potentials A a which are interpreted as the net momenta q of 
pairing in the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase. When l/v2 < Sfi/A < 1 where A is 
the gap energy, the effective potential suggests that the instability leads to a LOFF-like state which 
is characterized by color-rotated phase oscillations with small q. In the vicinity of 8fj,/A = l/y/2 
the magnitude of q continuously increases from zero as the effective potential has negative larger 
curvature at vanishing A a that is the Meissner mass squared. In the gapless 2SC (g2SC) phase, in 
contrast, the effective potential has a minimum at gA a ~ S/j, ~ A even when the negative Meissner 
mass squared is infinitesimally small. Our results imply that the chromomagnetic instability found 
in the gapless phase drives the system toward the LOFF state with q ~ Sfi. 

PACS numbers: 12.38.-t, 12.38.Aw 



Quark matter has a rich phase structure in the high 
baryon or quark density region. In a decade we have wit- 
nessed tremendous developments in theory, particularly 
in superconductivity of quark matter pj. Color super- 
conductivity is inevitable from the Cooper instability in 
cold and dense quark matter. In the asymptotic den- 
sity where the perturbative technique is applicable, the 
color-flavor locked (CFL) phase [2| where all quarks are 
gapped is concluded from the first-principle calculations 
of Quantum Chromodynamics (QCD). 

The lower density region we explore, the more compli- 
cated phase possibilities we have to encounter. The main 
reason why the situation is perplexing at intermediate 
density is that a "stress" between quarks which would 
form a Cooper pair is substantial when the quark chem- 
ical potential, \i q , is comparable to the strange quark 
mass, M s . Such an energy cost by the stress, or the Fermi 
energy mismatch 5[A, is necessary to bind two quarks into 
a pair with zero net momentum, q = 0. The stress can 
be reduced by making a pair between quarks sitting on 
different Fermi surfaces, which results in q ^ 0. If the 
energy gain by easing the stress is greater than the ki- 
netic energy loss coming from nonzero net momentum, 
the color superconducting phase with q ^ would be 
realized. Since such a state breaks rotational symmetry, 
this crystalline color superconducting phase |3(, that is, 
a QCD analogue of the Larkin-Ovchinnikov-Fulde-Ferrel 
(LOFF) phase Q , takes a crystal structure [j| . 

There is another different possibility to consider while 
keeping q — 0; once Sfx exceeds the gap energy, A, the 
Cooper pair tends to decay into two quarks. In other 
words the corresponding quarks have the energy disper- 
sion relation which is gapless. Such a phase is called the 
gapless superconducting phase |7]> 13 j that is, a QCD 
analogue of the Sarma phase 9] . It would need a careful 
comparison of energies to see which is favored in real- 
ity El- Interestingly enough, recently, these two 
different candidates, the crystalline and gapless super- 
conducting phases, have turned out to be closely related 




LOFF-like? Mixed Phase 




LOFF Mixed Phase 

FIG. 1: Schematic energy landscape: (a) The state falls from 
the unstable gapless phase down toward a LOFF-like phase 
which is separated from the LOFF phase (if it exists). There 
are alternatives such as the LOFF phase, mixed phase, and so 
on, which should be energetically compared to the LOFF-like 
phase, (b) The instability directly leads to the LOFF phase. 



through instability. 

The gapless superconducting phase is known to be un- 
stable in fact and has to give way to some other sta- 
ble states. In QCD the negative (color) Meissner mass 
squared exhibits what is commonl y re ferred to as the 
chromomagnetic instability [l2L fl3l HiL Il5| . This is the 
central issue we address. The chromomagnetic instability 
is to be interpreted as instability toward the single plane- 
wave LOFF state ^(J, as we will closely discuss later. 
Also it has been revealed that the two-flavor and three- 
flavor LOFF phases are chromomagnetic stable 0, 0] 
(see also Ref. However, it does not necessarily mean 

that the instability problem has already been resolved. 
The question we raise is as follows; can we simply iden- 
tify the instability-induced state with the stable LOFF 
phase? It is certain that the instability tends to favor a 
LOFF-like state, but such a LOFF-like state might exist 
separated from the LOFF phase (if it exists) as sketched 
in Fig. ma), an d tnen the proposed stable LOFF phase 
is not the destination from the instability but one al- 
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ternative free from the instability like a mixed phase. Of 
course, it might be possible that the instability is directly 
connected to the LOFF phase as sketched in Fig.QJb). 

To address this question we have to define the qualita- 
tive difference between the LOFF and LOFF-like phases. 
We shall distinguish them by their characteristic wave 
numbers. That is, if the net momentum q is given as of 
order Sfi, then we regard the system as going to the con- 
ventional LOFF state that has |q| ~ 1.25fi. If q is small 
enough to be well separated from q ~ Sfi inherent to 
the LOFF phase, we consider that the system is then in 
the LOFF-like state. For the purpose of clarifying which 
situation of Figs. [Ja) and (b) is more relevant, we will 
calculate the free energy as a function of q, or the color 
vector potential A Q (a being the adjoint color index) 
which is related to q through the covariant derivative. 

As a preparation for our discussions we shall briefly 
look over the color superconducting phases of our in- 
terest, the Meissner mass in respective phases, and as- 
sociated chromomagnetic instability. The predominant 
pairing is 
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with ip — Cip T and a, b and i, j being the color and 
flavor indices, respectively. Under this color-flavor locked 
ansatz, Ax ^ 0, A 2 ^ 0, and A 3 ^ defines the CFL 
phase, while the two-flavor superconducting (2SC) phase 
has A3 ^ and Ai = A2 = 0, that means only ru-gd 
and rd-gu quarks make a pair. The gapless 2SC and 
CFL phases (abbreviated as the g2SC and gCFL phases) 
occur when Sfi « fi e /2 > Sfif = A3 and Sfi ~ Mg/2fi q > 
6 /if = Ai, respectively, where fi e is the electron chemical 
potential. In the single plane-wave LOFF ansatz the gap 
parameters are augmented as 



LOFF 



exp[-2iq v ■ x] A n . 



(2) 



The Meissner mass is the screening mass for transverse 
gauge fields. The individual mass is a quantity dependent 
on the gauge choice; we can arbitrarily shuffle eight gluon 
fields A\ , . . . , Ag by a gauge rotation. It should be noted, 
however, that the choice of the diquark condensate 
specifies a gauge direction and then the Meissner mass 
is uniquely determined. The finite Meissner mass arises 
associated with spontaneous symmetry breaking and the 
Higgs- Anderson mechanism in superconductors |20| . 

In the 2SC phase Ai, A2, and A3 remain massless and 
the rest A4, . . . , As earn a finite Meissner mass. The elec- 
tromagnetic field A 7 has mixing with As which ends up 
with two eigen-fields Ag and Ay. The system has sym- 
metry among A4, . . . , A 7 and thus the resulting Meiss- 
ner mass is common for them, while Ag has a different 
Meissner mass. Because the modified (i.e. color-mixed) 
electromagnetic U(l) symmetry is unbroken, A 7 stays 
massless. It is known in the CFL phase at M s = 0, on 
the other hand, that all eight gluons Ai, . . . , A 8 have a 
common and nonvanishing Meissner mass. In the pres- 
ence of mixing with A 7 the massive eigen-field Ag pulls 
away from the others and A 7 is massless. 
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TABLE I: Chromomagnetic instability for each gluon in the 
two-flavor and three-flavor cases at zero temperature. Here 
Sfi^' 7 — Sfi^/y/2 and Sfii' 5 ^ Sfif' 7 has been numerically esti- 
mated as ~ 2.3 5/if in Ref. Jla| . 



In the two-flavor case at finite 5fi away from the ideal 
2SC phase, there are still only two independent Meissner 
masses |l2j; one is for A4, . . . , A7 and the other is for 
Ag. The Meissner mass squared for A4, . . . , Aj becomes 
negative (i.e. the Meissner mass is imaginary) for Sfi > 
Sfif 7 = SfiS/y/2, that means an instability occurs not 
only in the g2SC phase (Sfi > Sfif ) but in the 2SC phase 
(Sfi^T 7 < Sfi < <5/if) also. The Meissner mass squared for 
Ag is negatively divergent at the gapless onset, 6fi = 5fif, 
and remains on negative in the entire g2SC side. 

The three-flavor case with finite M s has a more com- 
plicated pattern and there are five independent Meissner 
masses hj. With nonzero M s one should take account 
of mixing among A3 , As , and A 7 properly, from which 
two massive eigen-fields, Ag and Ag, and one massless 
A 7 result. The Meissner mass is degenerated for A x and 
A2 due to symmetry, and so is for A4 and A5, and for 
Aq and A7. No instability takes place until the system 
reaches the gCFL phase. At Sfi — 5fi& negatively diver- 
gent Meissner masses squared appear for A1-A2 and for 
Ag. As for A4-A5 and A6-A7, when Sfi gets larger than 
critical values Sfi^~ 5 and 8fi®~ 7 respectively, they even- 
tually have negative Meissner masses squared, which is 
presumably related to the instability in the two-flavor 
calculation. We shall summarize the instability patterns 
in Table IJ 

Now let us consider what the negative Meissner mass 
squared signifies. It is a textbook knowledge that in 
the <p -theory, for the simplest example, a nonzero ex- 
pectation value (4>) 7^ grows when the screening mass 
squared for tp is negative. In the language of the effective 
potential the negative mass squared means that a state 
lies in a maximum of the potential and a true ground 
state should exist somewhere down away from (0) = 0. 
Therefore it is quite natural to anticipate that the chro- 
momagnetic instability is cured by nonzero color vector 
potentials (A Q ). Actually the Meissner mass squared is 
the coefficient of the quadratic terms in the kinetic energy 
expanded in A a , 
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where the kinetic energy term takes a form of 



rt kin [A,A] 

= K^,[(d% r ,, + igA™ l )A* v ,][(d% 
+ higher-order terms in A 



(4) 



due to symmetry. In the following discussions we shall 
ignore mixing with A 7 for simplicity. Here we should 
remark that the stiffness parameter k depends on the 
"flavor" indices rf and a' as well as the "color" indices r\ 
and a. This assignment is understood from that r\ of A v 
contains the information on both color and flavor as fixed 
in Q). The covariant derivative acting on a color-triplet 
rotates rj 1 into 77, while the flavor is intact as 77'. 

In the two-flavor case only K33 is relevant and we can 
forget about "flavor" , as parametrized in Refs. 0, |2^| . 
Then the chromomagnetic instability with two indepen- 
dent Meissner masses in this case can be expressed by 
a combination of two parameters k^ 1 ) and where 
K 33 = K ^$v° + K ^A^A* which makes fikin a color- 
singlet. The single plane- wave LOFF state character- 
ized by |J2j with only A3 nonvanishing is sensitive to 
+ k( 2 )|A3| 2 and thus such a gap parameter 
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ansatz cannot separate two distinct instabilities for As 
with m 2 MSS cx K33 and for A4, . . . , A7 with to|^ 4 . 7 cx 



if they coexist. 

It is intriguing to look into the three-flavor case next. 
The stiffness parameter can be decomposed as 



K r)'cT' — K os\ 6 m l^V^V'^' + K dia E 'W ^r^' 
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This decomposition is justified by the color and flavor 
structure in the quark one- loop calculations 0. Then 
the first term involving is relevant to the Meissner 
mass squared for the color-off-diagonal gluons; 



"4i- 2 ««off(|Ai| 2 + |A 2 | 2 ), 
^^offdAsP + IAil 2 ), 
m 2 Me . 7 « ^ ff (|A 2 | 2 + |A 3 | 2 ). 
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The Meissner mass squared for the diagonal gluons, on 
the other hand, comes from the second term involv- 

and m\ ISS are writ- 



ing K Lg' that is ' TO M33' 



'A/38' 



ten as a linear combination of six components of the 
symmetric 3x3 matrix fc^iag- We- note that, when the 
ansatz is substituted into Q, the instability toward 
finite q v does not reflect the information of k* s , that is, 
d 2 i} kin /dq r] dq <T oc «^f ag . 

From Q we can immediately understand how the chro- 
momagnetic instability is generally transformed to the 
instability toward the single plane-wave LOFF phase, 
which has been shown for As in explicit calculations in 
two-flavor quark matter in Rcf. 16] . The point is that 
the gauge fields in the quark sector appear only in the 
covariant derivative, so that they can be absorbed as a 
phase factor of the gap parameters. 

Now we assume that we have an instability only for Ag,. 
The rotational symmetry is broken and we choose the 



n-direction in three-dimensional spatial space in which 
Ag acquires an expectation value. Then the covariant 
derivative is equivalently rewritten as 



[8%^ - 5 in igA^(t 8 )^]A v > 



exp[igt 8 A$ ■ x] vv > d l {exp[-igt 8 A$ ■ x] v > v » A v „ }, 



(7) 



where rj a 's are the color group generators in the funda- 
mental representation. The color rotation results in 



/ exph-^Ag-x] AA 



exp[-igt 8 A 8 ■ x] ■ A = 
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exp 

y exp[+^|A 8 ■ x] A 3 J 



(8) 



which is nothing but the diquark condensate peculiar to 
the three-flavor single plane-wave LOFF state. [We as- 
sumed that A a is a constant, but the generalization to 
inhomogeneous A a (x) |23J is easy; the exponential part 
is then the Wilson line.] From the above rewriting, it is 
apparent that the non-LOFF (ordinary) superconduct- 
ing phase with a color vector potential Ag is equivalent 
to the LOFF phase whose spatial oscillation is charac- 
terized by A§ with no vector potential. Of course, this 
general argument works in the two-flavor case as well; 
Ai = A2 = and a phase factor emerges for A3 alone, 
so one could interpret such an overall phase as associated 
with the baryon number [24] , though such an interpreta- 
tion has only a limited meaning. 

One has to be careful when this argument is applied 
for the off-diagonal gluons, Ai, A 2 , A4, A5, Aq, and 
A7. For instance, if the instability occurs in A4, then 
the phase factor is no longer in the form of the single 
plane-wave. In the same way as in the previous case we 
have 



exp[— igt A A± ■ x] ■ A 

-i sin[| A 4 ■ x] A 3 + cos[| A 4 • x] Ai 
A 2 

cos[§ A 4 ■ x] A3 - i sin[§ A 4 • x] Ai 



(9) 



This represents not a single but rather multiple plane- 
wave LOFF state, or color-rotated single plane-wave 
LOFF. In the two-flavor case we keep A3 alone and then, 
interestingly, © indicates that we definitely need to have 
not only the third component cos[|A 4 -cc] A3 but also the 
first component — isin[|A 4 ■ x] A3 which is not consid- 
ered at all in the conventional two- flavor treatment. Our 
analysis agrees with the conclusion of Ref. that the 
single plane- wave LOFF state would still have instability 
for the off-diagonal gluons. 

From the discussions so far we can establish the qual- 
itative (apart from a color rotation) correspondence be- 
tween the vector potentials and the net momenta of the 
single plane- wave ansatz as 
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FIG. 2: Free energy difference from the energy without vector 
potentials in two-flavor quark matter as a function of A4. 
The critical Sfi for A4 is Sfi^T 7 = 5/if/V2 and the g2SC phase 
occurs at 5fi — 6fj%. 
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FIG. 3: Free energy difference from the energy without vector 
potentials in two-flavor quark matter as a function of As in 
the unit of A. 



Now we shall estimate the magnitude of characteristic q 
as a result of the instability using the above relations. 

For that purpose we need to know the higher-order 
terms in A a in the expansion J2J. As we will explain 
shortly, however, such an expansion is no longer valid 
in the gapless phase. Thus we must evaluate the A Q - 
dependent part of the free energy without expansion. For 
simplicity we will limit our discussions only to the two- 
flavor calculations from now on. 

We write down the 48 x 48 (two-flavors, three-colors, 
two-spins, particle-antiparticle, and two-Nambu-Gorkov- 
doublers) quasi-quark propagator with either A4 (that 




FIG. 4: Diagrammatic expansion of the free energy in terms 
of g{A a ) /e(p) where the quark energy e(p) comes from the 
quark propagator. 



we arbitrarily chose among A4, . . . Aj) or As and cal- 
culate the quasi-quark energy e(p) which depends on the 
momentum angle to the vector potential, i.e., p-A a . The 
free energy is available as integration of the sum over all 
48 |e(p)|'s with respect to p. We regulate the momentum 
integration by the ultraviolet cutoff A = 1 GeV and sub- 
tract the free energy at A = Sfi = to get rid of the cutoff 
artifact. The gap parameter is fixed at A = 100 MeV. It 
should be noted that the analytical formulae utilized in 
the two-flavor LOFF calculations do not work for 

A4 which is color off-diagonal. 

We present the numerical results in Figs. [5] and [3J The 
potential curvature at A a = corresponds to the Meiss- 
ner screening mass squared. In Fig. [3 the Meissner mass 
is real finite for Sfi < Sfi*~ 7 , while the origin A4 = be- 
comes unstable when Sfi > Sfi^.' 7 , as is manifest from the 
results at Sfi — Q.9Sfif 7 (dotted curve) and Sfi = l.lSfif 7 
(solid curve). This instability occurs continuously and 
we can see from the Sfi = l.2Sfif 7 results (dashed curve) 
that the expected A4 grows as Sfi approaches 5fif. It 
is known [l^ that the negative Meissner mass for A4 
becomes small again when Sfi is larger than 5/if. Cer- 
tainly our calculations for Sfi — 2.05 fif 7 = lASfif (dot- 
dashed curve) result in smaller potential curvature and 
thus smaller Meissner mass than those for Sfi = 1.2<5/i*" 7 . 
Nevertheless, the expected A4 is larger and we find a po- 
tential minimum at g|v44|/4 ~ 1.39A = 0.985/i. In this 
way the results for <5/i > <5/if make a sharp contrast to 
the nature of the instability for Sfi ~ S/i^' 7 < Sfxf. In 
the gapless region where Sfj. > 6 [if, the expected g|A 4 |/4 
(and thus q) is of order <5/i however small the potential 
curvature (Meissner mass squared) is. 

The same observation is apparent also in Fig. [21 The 
Meissner mass squared for As is negative divergent at 
5 fx = <5/if, meaning that the potential has a cusp at 
A% = then, which is confirmed in our results as seen 
at Sfi — 1.01<5/if (short-dashed curve). When Sfi pulls 
away from the onset value, the negative Meissner mass 
squared becomes smaller, and at the same time, Ag, ac- 
quires a larger expectation value of order Sfi again. 

We would emphasize that these findings are unex- 
pected results; if the Ginzburg-Landau expansion of Q 
works with a positive definite quartic term in A a , an in- 
finitesimal negative k (potential curvature) simply leads 
to an infinitesimal A a . Therefore, our results imply that 
not only the quadratic term but also quartic and even 
higher-order terms are significantly affected by gapless 
quarks when Sfi > 5/if. The reason why the Ginzburg- 
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Landau expansion breaks down can be understood in a 
diagrammatic way. 

Figure^Jshows an example of the diagrammatic expan- 
sion of the free energy in terms of A a . The dimension- 
less expansion parameter is obviously g(A a ) /e(p) where 
e(p) is the quark energy stemming from the propagator. 
Therefore such an expansion is no longer legitimate once 
gapless quarks whose e(p) can become vanishingly small 
enter the loop. In other words, in the gapless phase, the 
Meissner mass squared is far from informative on the true 
ground state; the smallness of the Meissner mass squared 
does not mean the weakness of the instability. 

It should be noted that our potential analysis neither 
preserves neutrality nor solves the equations of motion. 
As long as A a is small, the coupling between A a and the 
other parameters such as the chemical potentials and the 
gap parameters is small due to approximate rotational 
symmetry. We can thus expect that the free energy we 
estimated would not be modified significantly by neutral- 
ity and the condensation energy in the region where A a 
is small. This is not the case, however, once A a ~ <5/j. de- 
velops. Hence, strictly speaking, we cannot say anything 
about the exact location of the potential minimum, but 
we can at least insist that there is no stable LOFF-like 
state with small q directly resulting from the chromo- 
magnetic instability in the gapless phase even when the 
negative Meissner mass squared is tiny. 

Our results in the two-flavor case are not direct evi- 
dence but suggestive; we can anticipate that the situa- 
tion in Fig. ^a) is realized for the off-diagonal gluons 
at 5 He' 7 < 6 (j. <C 6[j,%, while the situation in Fig. IHb) is 
likely to be the case for 6/j, > <5/x§. In the three-flavor 
case we conjecture that the instability picture is close 
to Fig. [Hb) since the chromomagnetic instability then 
occurs only in the gapless region of Sfi. These are our 
central conclusions derived from the numerical results. 

Finally let us comment on the possibility of coexistence 
of both A4 and Ag in the two-flavor case. We shall call 
such a state the gluonic phase |25|. One should be careful 
about the terminology not to fall in a mere interpreta- 
tion; we would use the nomenclature, the gluonic phase, 
differently from the original usage in Ref. [2j|, but to 
mean a state in which all A a in the covariant derivative 
cannot be simultaneously removed by any gauge rotation 
of A^. Thus one can uniquely define the gluonic phase 
in a way distinct from the LOFF and LOFF-like states. 

If A4 and Ag are not parallel, A\ and A\ for in- 
stance, then we cannot find an appropriate gauge rota- 
tion A — * VA to eliminate them simultaneously. Namely, 



the gauge rotation matrix V satisfying V'd 1 V — igA\t A 
and V^d 2 V = igAgt 8 does not exist if V is assumed not 
to have any singularity. 

The gluonic phase has one more significant difference 
from the LOFF state besides the covariant derivative; it 
has nonvanishing chromomagnetic field. For our example 
A\ and A\ produce a nonzero field strength tensor, 

Bl = F al2 = -^-5 a5 gA\Al (11) 
which means that the system has a uniform chromomag- 
netic field in it. However, such a state would never be 
realized, otherwise the field energy diverges. To put it in 
another way, the vector potentials A\ and A\ do not solve 
the Yang-Mills equations of motion, D^F^ a = |2(j|. 
Therefore, we do not think that the gluonic phase results 
from the chromomagnetic instability. 

In fact, one can reduce the field energy by making A4 
and Ag be parallel to each other. Then the Yang-Mills 
action simply vanishes. This argument can be easily ex- 
tended to more generic A a in the three-flavor case. We 
would thus reach a conclusion that all nonvanishing A a 
as a result of the chromomagnetic instability are aligned 
to the same direction energetically. Then such A a can be 
eliminated by a gauge rotation of A,,. That is, the likely 
destination is a state characterized by the gap parame- 
ters exp[— igt a A a ■ x] • A where the summation over a is 
taken. This is what is called the colored crystalline phase 
in Ref. and, as we have seen in ©, characterized 

by a LOFF ansatz beyond the single plane-wave one. Al- 
though the difference from the single plane-wave ansatz 
exp[— 2iq v ■ x] A v is just a color rotation, it changes the 
physics because rj of A v has the information of "flavor" 
as we have already discussed in kI^i ■ 

In summary, based on our numerical results for the 
free energy as a function of the color vector potentials, 
we have reached a speculation that the chromomagnetic 
instability in the gapless color superconducting region 
leads to the LOFF phase, meaning that the LOFF phase 
is not an alternative but a destination of the instability. 
In contrast, the instability found in the 2SC (not g2SC) 
phase drives the system toward a LOFF-like state which 
is qualitatively distinct from the LOFF phase. 
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